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One of the most exciting features of the recently dis¬ 
covered double giant dipole resonance (DGDR) is the ab¬ 
solute value of its excitation cross section in relativistic 
heavy ion collisions. These values were extracted from 
total cross sections by separating a contribution arising 
from the excitation of single giant dipole (GDR) and 
quadrupole (GQR) resonances and were found to be en¬ 
hanced by factors of order of 2-3 for ^^®Xe and ^®^Au 
1^ as compared to any theoretical calculations available. 
Although for ^°®Pb the experiment-theory correspon¬ 
dence is much better, theoretical calculations still under¬ 
predict the DGDR cross section by about 30%. 

At the present time we have an experimental proof 
1^ that Goulomb excitation of the DGDR in relativis¬ 
tic heavy ion collisions occurs in a two-step process. The 
DGDR states are embedded in a sea of other two-phonon 
states. In this paper we consider the sea contribution to 
the total cross sections. It will be concluded that numer¬ 
ous two-phonon states, other than [GDR G GDR] ones, 
excited in one-step processes form a physical background 
in the DGDR energy region which has to be taken into 
account in the analysis of experimental data. 

The first investigation of the direct excitation of two- 
phonon states at high energies from the ground state was 
done in refs. |^,^. The main purpose of these papers was 
to look for resonance structures in photoexcitation cross 
sections related to the excitation of rather selected two- 
phonon configurations. Although in ref. it was pointed 
out that [GDR 0 GQRji- and [GDR G 2]^]i- configura¬ 
tions were not the only ones (among many different two- 
phonon [1“ G 2+]i- states) to determine the total cross 
section, two-phonon 1“ configurations made of phonons 
of other multipolarities were omitted in calculation. 

While dealing with electromagnetic, or with Goulomb 
excitation from a O'*" ground state, the priority attention 
has to be paid to the final states with the total angu¬ 
lar momentum and parity = 1~. Making use of the 
formalism presented in ref. |^] we have calculated firstly 
the cross section of photoexcitation of two-phonon states 
[A^^ G where A^^ and AJ^ are both natural A’^" 

(tt" = (—I)'*') and unnatural A'^’* (7r“ = (—1)^+^) par¬ 
ity phonons with multipolarity A from 0 to 9. A phonon 
basis is obtained by solving quasiparticle-RPA equations 
for each multipolarity A’^ within the quasiparticle-phonon 


model (QPM) [|j. These equations provide a set of 
one-phonon states X^(i) with the same spin and parity, 
but with different excitation energies E(i) and internal 
fermion structure of phonons; the index i is introduced 
to distinguish between them. 

The results of the calculation for ^^®Xe and ^°®Pb in¬ 
tegrated over the energy interval from 20 to 35 MeV are 
presented in table |. Each configuration [A^^ G A^^] in 
the table means a sum over a plenty of two-phonon states 
made of phonons with a given spin and parity A]^^, A^^, 
but different RPA root numbers fi, i 2 of its constituents 

ct([Ai^ g A2"]) = ® • 

il,i2 

The total number of two-phonon 1“ states included in 
this calculation for each nucleus is about 10^. The ab¬ 
solute value of the photoexcitation of any two-phonon 
state under consideration is negligibly small but alto¬ 
gether they produce a sizable cross section. Table | 
demonstrates that different two-phonon configurations 
give comparable contributions to the total cross section 
which decreases only for very high spins because of the 
lower densities of such states. As a rule, unnatural par¬ 
ity phonons play less important role than natural parity 
ones. For these reasons we presented in the table only the 
sums for [naturalGunnatural] and [unnaturalGunnatural] 
two-phonon configurations. 

The cross section of the photoexcitation of all two- 
phonon 1“ states in the energy region 20-35 MeV from 
the ground state equals in our calculation to 511 mb 
and 423 mb for ^^®Xe and ^°®Pb, respectively. It is 
not surprising that we got a larger value for ^^®Xe than 
for ^°®Pb. This is because the phonon states in Xe are 
composed of a larger number of two-quasiparticle con¬ 
figurations due to the pairing. The same values for 
two-phonon states with angular momentum and parity 
J'^ = 2+ are an order of magnitude smaller. We point 
out that the direct excitation of [l“Gil“] 2 + configuration 
or [GDR ® GDR] 2 + is negligibly weak. The calculated 
values should be compared to the cross section of the 
photoexcitation of the single-phonon GDR which in our 
calculation equals to 2006 mb and 2790 mb, respectively. 
A contribution of two-phonon 1“ states to the total cross 
section at GDR energies is weaker than at higher ener¬ 
gies because of the lower density of two-phonon states 
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and lower excitation energy and can be neglected consid¬ 
ering the GDR itself. 

For photoexcitation cross sections are known 

from experimental studies in ( 7 , n) reactions up to the 
excitation energy about 25 MeV ||^^. It was shown that 
QPM provides a very good description of the experimen¬ 
tal data in the GDR region 0, while theoretical calcu¬ 
lations at higher excitation energies which account for 
contributions from the single-phonon GDR and GQR™ 
essentially underestimated the experimental cross section 
1^. The experimental cross sections above 17 MeV are 
shown in Fig. together with theoretical predictions. 
The results of the calculations are presented as strength 
functions obtained with averaging parameter equal to 
1 MeV. The contribution to the total cross section of 
the GQRi„ (short-dashed curve), the high energy tail 
of GDR (long-dashed curve), and their sum (squared 
curve), are taken from ref. |^. The curve with trian¬ 
gles represents the contribution of the direct excitation 
of the two-phonon states from our present studies. The 
two-phonon states form practically a flat background in 
the whole energy region under consideration. Summing 
together the photoexcitation cross sections of all one- and 
two-phonon states we get a solid curve which is in a very 
good agreement with the experimental data. 

From our investigation of photoexcitation cross sec¬ 
tions we conclude that in this reaction very many dif¬ 
ferent two-phonon states above the GDR contribute on 
a comparable level, forming altogether a flat physical 
background which should be taken into account in the 
description of experimental data. On the other hand, 
Coulomb excitation in relativistic heavy ion collisions 
provides a unique opportunity to excite a very selected 
number of two-phonon states by the absorption of two 
virtual 7 ’s in a single process of projectile-target interac¬ 
tion 1^. Theoretically this process is described using a 
the second order perturbation theory of the semi-classical 
approach of A. Winther and K. Alder [^0|j|]. Since ex¬ 
citation cross sections to second order are much weaker 
than to first order of the theory, only two-phonon states 
connected to the ground states by two El-transitions can 
be observed. These two-phonon states have the structure 
[l“(i) G l~(i 0 ]o+. 2 + and form the DGDR. 

To describe the properties of the DGDR in ^°®Pb we 


applied the technique developed in refs. 00 - Within 
this technique we couple two-phonon [l“(*)( 8 il“(i 0 ]o+, 2 + 
states with many three-phonon ones. Two-phonon states 
are built from the same six l“(i) RPA phonons as in 
ref. [|^ which have the largest B(E1) values and exhaust 
90.6% of the EWSR. Only 0"^ and 2+ components of the 
DGDR were considered (see ref. for the quenching 
of the I'*' component). We included in the calculation 
about 7000 three-phonon states which have the largest 
matrix elements for the interaction with the selected 21 
two-phonon states. Diagonalization of the QPM Hamil¬ 
tonian on the basis of these two- and three-phonon states 
yields a set of 0 + (and 2 +) states; the wave function of 
each state includes all two- and three-phonon terms with 
different weights for different states. To distinguish be¬ 
tween these states we introduce the index v. To obtain 
the Goulomb excitation cross section in second order per¬ 
turbation theory we also need to have the structure of the 
GDR as an intermediate state. For that we calculated the 
GDR fine structure by coupling the same six strongest 
one-phonon l“(i) states to about 1200 two-phonon 1 “ 
states in the GDR region. 

The cross section of the DGDR(i/o+, 2 +) states excita¬ 
tion via the GDR(i/i-) states in this reaction equals to 

“ I ,+ ) < 1 (^1- )II'^1||0^S. ^ 

X < [1-®1-](j.o+.2+)I|A1||1-(:/i-) > I' 

where A{Ei,E 2 ) is the energy dependent reaction am¬ 
plitude. A straightforward calculation of the two-step 
process for the excitation of 7000 DGDR states via 1200 
intermediate GDR states is very time consuming. Mak¬ 
ing use of a very smooth dependence of the function 
A{Ei,E 2 ) on both arguments we tabulated this func¬ 
tion and used it in the final calculation of the DGDR 
Coulomb excitation cross section in relativistic heavy ion 
collisions. We considered the excitation of the DGDR 
in the projectile for a ^°®Pb (640 A-MeV) -|-^°®Pb colli¬ 
sion, according to the experiment in ref. and used the 
minimum value of the impact parameter, b= 15.54 fm, 
corresponding to the parameterization of ref. . 

The cross section for Coulomb excitation of the DGDR 
is presented in Fig. || by the short-dashed curve as a 
strength function calculated with an averaging parame¬ 
ter equal to 1 MeV. The width of the DGDR for ^°®Pb is 
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very close to \f2 times the width of the single GDR, as 
for the case of ^^®Xe [^. The contribution of the back¬ 
ground of the two-phonon 1 “ states to the total cross 
section is shown by a long-dashed curve in the same fig¬ 
ure. It was calculated in hrst order perturbation theory. 
The role of the background in this reaction is much less 
important than in photoexcitation studies. First, it is 
because in heavy ion collisions we have a special mech¬ 
anism to excite selected two-phonon states in the two- 
step process. Second, the Coulomb excitation amplitude 
is exponentially decreasing with the excitation energy, 
while the El photoexcitation amplitude is linearly in¬ 
creasing. Nonetheless, Fig. ^ shows that the direct ex¬ 
citation of two-phonon 1 “ states cannot be completely 
excluded from consideration of this reaction. Integrated 
over the energy interval from 20 to 35 MeV these states 
give a cross section of 50.3 mb which should be compared 
to the experimental cross section in the DGDR region for 
the ^°®Pb (640 A-MeV) -|-^°®Pb reaction which is equal 
to 380 mb 1^. Taking into account that the reported 
experiment/theory enhancement is of about 30% for this 
reaction, the two-phonon 1 “ states omitted in extract¬ 
ing the DGDR strength from experimental cross section 
maybe responsible for an appreciable part of the effect. 
Appreciable values of one-step processes in DGDR exci¬ 
tation is not in contradiction with experimental findings. 
It is known that Coulomb excitation of a projectile in an 
n-step process has the following dependence on the target 
charge: The reported value n = 1.8(3) in the 

DGDR region for ^°®Pb allows for some contribution 
of one-step transitions. 

The solid line in Fig. ^ is the sum of DGDR and two- 
phonon background excitations in relativistic heavy ion 
collisions. As mentioned above, this curve has a closer 
correspondence to the experimental cross sections re¬ 
ported for the DGDR than the short-dashed one. Cen¬ 
troids and widths calculated for the short-dashed and 
solid curves of this figure display other interesting fea¬ 
tures. Taking into account the background of two- 
phonon 1“ states results in a visual shift of the DGDR 
centroid by -200 keV and a 16% increase of the “DGDR” 
width. The same effect, although with large experimental 
uncertainties, was reported in an experiment with respect 
to the harmonic picture of the DGDR excitation j|]. We 


point out that the effect of the visual shift of the DGDR 
centroid is even a bit larger than anharmonicity effects 
examined for this nucleus in ref. [T^ . 

The effect of the direct excitation of two-phonon states 
from the ground state was considered in the present paper 
under the assumption that these states do not interact 
between each other and with other one- and three-phonon 
configurations. The first type of interaction, unharmonic- 
ity effect, is rather weak as discussed in the refs, 15 


and cannot change our conclusions as far as the wide en¬ 
ergy region is considered. But one may argue that a weak 
coupling to one-phonon 1 “ states with much larger ma¬ 
trix elements of photoexcitation from the ground state 
may reduce the calculated cross sections due to the de¬ 
structive interference with our two-phonon 1 “ states. Of 
course, the total contribution of this interference term 
integrated from 0 MeV to infinity equals to zero because 
of orthogonality of wave functions but it cannot be ex¬ 
cluded from a consideration for a limited energy inter¬ 
val. To estimate its influence on our results we have per¬ 
formed a calculation in which 1 “ states in ^°®Pb where 
described by the wave function consisting one- and two- 
phonon components. The interaction between these com¬ 
ponents has been calculated microscopically ^ . The full 
calculation with including all 10 ^ two-phonon 1 “ config¬ 
urations is not possible and we have done several tests 
with different truncations of the two-phonon basis ending 
up in total 10^ components for each test. All one-phonon 
configurations with the energy below 35 MeV have been 
included in calculations. Selecting important two-phonon 
configurations we have chosen the ones which have the 
largest matrix elements of the electromagnetic excitation 
< [Aj^(ii)( 8 )A 2 ^(z 2 )]i-||El||g.s. > and the largest matrix 
elements of interaction with one-phonon 1“ states. The 
results of our test calculations are the following. 

We do observe the effect of interference between tran¬ 
sitions to one- and two-phonon components of the wave 
functions of 1“ states. In all runs the interference is 
constructive in the whole energy region under consider¬ 
ation leading to an additional enhancement of the cal¬ 
culated cross sections. It becomes destructive only be¬ 
low 14.5 MeV reducing a little the total El-strength 
for the GDR. In the DGDR energy region we estimate 
from our calculations the size of the interference effect 
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roughly as 30-50% of what we get from the direct ex¬ 
citation of non-interacting two-phonon 1“ states. We 
think that this estimate is an upper limit for be¬ 

cause many other two-phonon configurations will block 
somehow a coupling of selected two-phonon configura¬ 
tions to one-phonon ones. On the other hand, the in¬ 
terference between transitions to one- and two-phonon 
configurations may play much important role for ^^®Xe 
in which the interaction between these configurations is 
sufficiently stronger as compared to ^°®Pb. It requires an 
extra study. 

In any case, the influence of a sea of two-phonon 1“ 
states above 20 MeV on calculated cross sections con¬ 
sidered in the present paper in the approximation of the 
pure two-phonon configurations can be taken as a lower 
limit of the real effect. 

In conclusion, we investigated the contribution of the 
direct excitation of a sea of two-phonon states above 
the GDR to photo-neutron cross sections and to cross 
sections in relativistic heavy ion Coulomb excitation. 
These states form a flat structureless background which 
is very important for the correct description of the photo¬ 
neutron data. Due to the peculiarities of Coulomb exci¬ 
tation in heavy ion collisions, its role is less important in 
this case, but its consideration appreciably removes the 
disagreement between experiment and theory. 

We thank Dr. V.V. Voronov for his suggestion to 
check the role of the interference effect discussed at 
the end of the paper. This work was partially sup¬ 
ported by the RFBR (grant 95-02-05701) and the Brazil¬ 
ian funding agencies CAPES, CNPq, FUJB/UFRJ, and 
by MCT/FINEP/CNPQ(PRONEX) under contract No. 
41.96.0886.00. 
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TABLE 1. Cross sections for the direct photoexcitation of 
different two-phonon configurations from the ground state in¬ 
tegrated over the energy interval from 20 to 35 MeV in ^®®Xe 
and ^°®Pb. The GDR cross section integrated over the energy 
of its location is presented in the last line for a comparison. 



(T.y, mb 


Conhguration 

i36Xe 

208 pb 

[0+ o i-]i- 

4.4 

3.9 

[1- (g) 2+]i~ 

36.6 

44.8 

[2+ (g)3-]i.^ 

82.8 

33.1 

[3- ® 4+]i- 

101.0 

56.7 

[4+ (g)5-]i- 

68.9 

37.3 

[5- ® 6+]i- 

49.2 

46.2 

[6+ (g)7-]i~ 

31.9 

49.8 

[7- ® 8+]i- 

13.6 

12.5 

[8+ (g)9-]i~ 

4.9 

9.0 

E [Af ®Af]i- 

Ai, A2 = 1 

71.4 

58.5 

E [Af®A;=]i- 

Ai, A2=1 

46.7 

71.1 

E [A:"'“®Af\- 

Ai,A2=0 

511.4 

422.9 

[GDR® GDR]2+ 

0.33 

0.22 

E [a:"®A2^"]2+ 

Ai, A2 = 1 

38.1 

21.7 
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FIG. 1. Photo-neutron cross section for ^°®Pb. Experimen¬ 
tal data (dots with experimental errors) are from ref. . The 
long-dashed curve is the high energy tail of the GDR, the 
short-dashed curve is the GQRi„ and the curve with squares 
is their sum. The contribution of two-phonon states is plot¬ 
ted by a curve with triangles. The solid curve is the total 
calculated cross section. 

FIG. 2. The contribution for the excitation of two-phonon 
1“ states (long-dashed curve) in first order perturbation the¬ 
ory, and for two-phonon 0"^ and 2^ DGDR states in second 
order (short-dashed curve). The total cross section (for ^°®Pb 
(640 A-MeV) -|-^°®Pb) is shown by the solid curve. 
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